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Abstract 



An extended Anderson model, including screening channels (non- 
hybridizing, but interacting with the local orbit), is studied within the 
Anderson- Yuval approach, originally devised for the single-channel Kondo 
problem. By comparing the perturbation expansions of this model and a 
generalized resonant level model, the spin-spin correlation functions are cal- 
culated which show non-Fermi liquid exponent depending on the strength of 
the scattering potential. The relevance of this result to experiments in some 
heavy fermion systems is briefly discussed. 
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I. INTRODUCTION 



The Anderson model, originally proposed to explain the formation of the local magnetic 
moment in metals, has played a fundamental role in exploring the correlation effects in 
many-body systems. Various theoretical approaches, including the Hartree-Fock |IJ , scaling 
theory 0, numerical renormalization group [H, perturbation expansion Qj, variational cal- 
culations ||, slave-boson mean field theory II, large-N diagrammatic expansion |/J, exact 
Bethe-Ansatz solution ||, and many others [fj, have been applied to study this important 
issue. The general understanding has been that the Anderson model, in both electron-hole 
symmetric and asymmetric cases (including the valence fluctuation regime) exhibits a local 
Fermi liquid (FL) behavior, i.e., the impurity contribution to the specific heat and magnetic 
susceptibility is regular. 

Recently, this "common understanding" has been questioned by studies of a generalized 
Anderson model, including the so-called screening channels [fL0|-|13|l . In fact, it was realized 
long time ago by Anderson and Haldane that the original Anderson model is not complete to 
fulfil the local charge neutrality Friedel sum rule in the mixed- valence regime, i.e., to satisfy 
two different equations, corresponding to two different valence states by varying only one 



single parameter-the local state energy level jnj. In addition to the channel of conduction 
electrons which hybridize with the local state, "screening" channels have been included, 
which do not hybridize with the local state but are related to it via Coulomb interactions. 



Haldane has also presented a Hartree-Fock mean field theory for this extended model |14 
which, unfortunately, missed the non-FL behavior contained in it. 

Recently, Varma and collaborators [|H],[iT[] have revived interest to this model, considering 
it as a single site version of the three-band Hubbard-type model proposed to describe the 



cuprates fL5| . The "screening" channels correspond to oxygen orbits, not mixing with copper 
due to symmetry. Their original motivation was to find a microscopic justification for 
the phenomenological marginal FL theory [plj. At first, some numerical evidence of non- 
FL behavior was provided by a Wilson renormalization group study [110]. Later, a strong 



2 



coupling limit Hamiltonian was derived which exhibits a quantum critical point separating 



the Kondo and empty-orbit regions [|TT|,0. The correlation functions around this critical 
point show non-FL behavior. However, the calculation of the correlation functions is very 
delicate due to the presence of the single occupancy constraint. We believe the correct result 
has been obtained in [JHJ which demonstrates a power law singularity for the impurity specific 
heat and magnetic susceptibility, namely, Ci mp /Co ~ Ximp/xo ~ T -3 / 4 in the unitarity 



limit (phase shift = 7r/2), in contrast to the logarithmic singularity anticipated earlier [11 
Whether this result is relevant for the high T c materials, is an open question. However, 
some heavy fermion compounds show such type of singularities. In particular, the alloys 
UPd x Cu 5 _ x (x = 1, 1.5) have been studied in detail, using both static and neutron scattering 
techniques [0. The critical exponent extracted from the experimental data []T8| , |19l is A = 



1/3, which is qualitatively consistent but quantitatively different from the theoretical result 
obtained for the generalized Anderson model ]13| . 



The calculations of |13| have been carried out using the bosonization technique and 
the canonical transformations which provide values for the critical exponents only in the 
unitarity limit. In view of the importance of this issue and the need to obtain results away 



from the unitary limit |2(j[ it is highly desirable to have another way to reconfirm and extend 



the previous result | 13| . In this paper we apply the Anderson- Yuval (AY) approach [^TJ to 



consider the extended Anderson model. 

The AY approach [^TJ was originally devised to study the single channel Kondo problem. 
The main idea of this approach is based on the time-dependent one-body formulation of 
the X-ray edge problem, developed by Nozieres and De Dominicis (ND) [p2]| . Using the 
AY perturbation expansion Toulouse could map the Kondo problem onto a resonant level 



model from which he could derive the well-known strong-coupling Toulouse limit ||23|| . Later, 
the bosonization technique was used by Schlottmann [^4]] to calculate physical properties 
in this strong coupling limit, although bosonization was employed much earlier by Schotte 
to derive the AY perturbation expansion for the partition function P5|] . The equivalence 
of this perturbation expansion with that of a resonant level model (at and away from the 



Toulouse limit) was shown explicitly and made use of for studying the physical properties 
by Wiegmann and Finkel'stein |[26|| . Recently, Fabrizio, Gogolin and Nozieres have gen- 
eralized the AY approach to consider the asymmetric two-channel Kondo problem and the 
FL-non-FL crossover within that model. They have mapped term-by-term the perturbation 
expansion of the two-channel Kondo problem to that of a generalized resonant level model. 
Under the scaling assumption these authors could provide an analytical description of the 
FL-non-FL crossover. In this paper we will follow their approach rather closely. Bosoniza- 
tion and AY approach are similar in the sense of scaling assumption, but the latter is not 
limited to the Born approximation and can be used away from the unitarity limit. More- 
over, it is a "brute force" partial resummation of diagrams without resorting to canonical 
transformations used in the bosonization approach. 

We should mention that there is another generalization of the Anderson model considered 
by Si and Kotliar p8|j29[| , who have included all possible density-density interactions of the 



hybridizing channel without invoking the screening channels. Using the renormalization 
group expansion they considered the weak coupling case PSL while the strong coupling 



limit was treated by bosonization P5| . In the final section we will discuss how their results 
are related to ours in the overlapping region. 

The rest of the paper is organized as follows: The model is defined and the basic strategy 
to treat the model is given in Sec. II, while the perturbation expansion and the mapping 
onto a generalized resonant- level model is presented in Sec. III. Furthermore, the correlation 
functions are calculated in Sec. IV to be confronted with experimental results. Finally, some 
concluding remarks are made in Sec.V. 

II. THE MODEL AND THE BASIC STRATEGY 

The Hamiltonian we consider in this paper is given as a sum 

H T = H + H, 

where the hybridizing part is 



H = H + Hj + t$>+(OK + ^c CT (0)), 

a 

Hp = e k c^ a c k(7 , 



a,k 



Hi = v J24(o)c a (o)( n(T -h, (i) 

cr Z 

while the screening part is 

= £ eikcLcjfe + £ V^ C +(0)Q CT (0)(n - I). (2) 

l,k,cr l,a ~~ 

Here c^ CT , Ck a are conduction electron operators in the hybridizing channel, with c+(0), c a (0) 
as their Fourier transforms at the origin of the coordinates (where the impurity is lo- 
cated). Similarly, cf ka , cik a are conduction electron operators in the screening channels. 
n — J2a n (T, n a = d+d a is the local impurity operator. The Hubbard U on impurity itself 
has been taken as infinity, so the only allowed states are |0>, |cr >, cr =f, j. Vo, V\ are 
the Coulomb interactions of the local electron with conduction electrons in the hybridizing 
and screening channels, respectively. The essential part of this Hamiltonian is the same as 
that in earlier papers [|TT|-|i"3|]. We have not included here the anti-parallel spin Coulomb 
interactions Vq J2a c t c a( n a ~ |) an d the spin-flip scattering V± J2a c t c °dtd a in the hybridiz- 
ing channel, considered in P^|T3| ]. As shown in those references, these terms do not affect 
essentially the behavior in the strong coupling limit. For simplicity, we have also taken the 
local level at the Fermi level, since we are mainly interested in the critical behavior itself, 
rather than the level renormalization per se. 



Before proceeding, let us briefly recall the basic strategy of ND |22| and AY |2T| to see why 
their approach can be applied to our problem. The crucial term in the X-ray edge problem 
is V J2 c + (0)c(0)dd + , where d is the deep-level electron, while c + (0), c(0) are the conduction 
electron operators at the origin. As it stands, this is a many-body problem. However, ND 
have realized that it can be converted into a time-dependent one-body problem, because the 
scattering potential is effective only after X-ray absorption or before X-ray emission (when 
dd + = 1). Since the internal degrees of freedom for the local electron are not involved, 
its propagator in the standard many-body technique can be traced out, thus converting it 



into a one-body problem. Moreover, instead of the usual Fourier representation of Green's 
functions, ND solved the Dyson equation directly in the time domain. The integral equation 
then obtained turned out to be singular, of the Muskhelishvili-type [^j from which the 
power-law time asymptotic behavior is extracted with exponents expressed in terms of a 
phase shift S = tan -1 (70/9 V), where u is the conduction electron density of states. 

Soon afterwards AY [|2l| realized that the ND trick can be used to obtain the perturbation 



expansion for the Kondo problem with a Hamiltonian 

Hk = \jY,\ c t c l S ~ + 4 c t 5+ + (4 C T - c t c S ^ ( 3 ) 
where, S + , S~ , S z are local spin components, while ct, ... are conduction electron 
operators at the origin. For any given sequence of spin flips at moments t\, t n (from up 
to down), t' n (from down to up) in a perturbation expansion of the partition function, 

the asymptotic expression can be used if the time difference ti — t i7 is much greater than 
the transient time to ~ the inverse of the bandwidth. The exponents can again be expressed 
in terms of the phase shift S = 5 + — 5- = 2 tan" 1 2L |^, due to the difference of scattering 
potential experienced by the up and down spins. By mapping this perturbation onto an 
expansion for some kind of resonant level model which can be solved exactly in limiting 



cases, the low energy physics can be extracted. This was the basic strategy of |26j and f27 
and will be followed in this paper. Instead of | |> and | |> two states for the single channel 
Kondo problem, we have here three states |0 >, | |>, and | |>. Also, we have screening 
channels in addition to the hybridizing channels. Nevertheless, the above programme can 
be still implemented, as seen from the next Section. 

III. PERTURBATION EXPANSION AND MAPPING ONTO A GENERALIZED 

RESONANT LEVEL MODEL 



Like in the X-ray edge and Kondo problems, we are interested in the time evolution of 
the system described by 

F(t) =< 0\e iHTT \0 >, 
6 



where Ht is the total Hamiltonian given by ([I]) and (0), while |0 > is the ground state which 
is degenerate in the mixed valence regime with the local electron in one of the states \oc >, 
a = 0, | and j. We consider the perturbation expansion in terms of the hybridization 
parameter t 

F(t) =< 0\e iH ' T T{e l Io^' H ^^}\0 >= f] f ' dr 2n P dr 2n ^ ■■■ H dn 

n=0 Jo Jo Jo 

< 0\e iH ' {T - T2n \H h ■ ■ ■ e iH ' {T ^ +1 - T ^iH h e iH ' (T2 ^' l] iH h ■ ■ ■ iH h e iH ' Tl \0 > (4) 

with 

a 

H' = H + Hj + H. (5) 

Due to the presence of <i+, d a in Hh, only even order terms are kept. A typical term will 
contain either spin up (down) operators only, or mixed. However, in view of the single 
occupancy constraint the spin up and down states can be connected only via the empty 
state. Therefore, these terms can be separated into up and down blocks. Using the relations 

e iH cJi\ a > if U >= |o > 

d+e lHT *\a>={ (6) 
0, if | a >^ |0 >, 

■ 0, if \a \a >, 

d a e* H ^\a>={ , ^' (7) 

e lH ° T i\0 >, if \a >= \a > 

with H a =< >, we can trace out d a operator. Since the Hamiltonian H a is different 

for the neighboring time moments Tj and r J+ i (or Tj_i), the calculation is similar to the 
X-ray edge or Kondo problem and we can apply AY approach to handle it. Explicitly, we fix 
T\j (r-[ ■) as moments when a local spin- up electron is created (annihilated). Likewise r 2 j (r 2 „-) 
for local spin-down electron. As said before, due to the single occupancy constraint, the only 
allowed sequences are: (ry, r^), (r 2 j, r^), i.e., the up and down sequences do not intersect 
each other. 



In the absence of the scattering potentials Vq = V\ = 0, the propagator for the conduction 
electron in the local presentation G a {r) =< Tc a (r)c£(0) > is [ p2[ : 



G<°>(r) 



IISq 



(8) 



t ~ i£o ^gn-r' 

where £o is the bandwidth, serving as a cut-off. The contribution of 2n vertices connecting 



r al , r aj (a =1,2) will be |21|,^7| 



D„ 



Y[(r ai - T aj ) JJ(r^ - T aj ) 

i<3 i<3 



n< 



r„i - r. 



(9) 



In the absence of the scattering Vq = VJ = 0, the contribution of separate terms in (|J) 
will be 

U = D^D^ (10) 

Now we turn on the scattering potential. First consider the hybridizing channel. For the 
spin up conduction electron at moment tu a scattering potential is switched on to give a 
phase shift 



5 = 25 = 2tan- 1 (^; 



(11) 



while at moment r' u an opposite phase shift is produced. On the other hand, the spin down 
conduction electron does not experience any phase shift at these moments. Of course, at 
moments r 2 j and r' 2i the situation is reversed. There are two types of contributions to the 
renormalization of the conduction electron propagation: from an open line Ul and closed 
loops (vacuum fluctuations). As shown earlier plp7| , 



U L = {D^D{)- 2& I\ 



(12) 



In deriving ( 12 ) we have taken into account the fact that in the process |0 >^ | f>, nothing 
is changed for H a , with a = 0, j, so there are no crossing terms. 



S 



Next we consider the screening channels. For channel I the conduction electron gets a 
phase shift 

^ = 25(0) / = 2tan- 1 (^) (13) 

at both T\i and r 2 j, and an opposite phase shift at t u and r 2i . The screening electrons con- 
tribute only to the closed loops. However, as seen from (0), the screening channel electrons 
scatter on the total charge n = J2a n cn so the process | |>— > |0 > will affect also spin down 
conduction electron in the screening channels. Therefore we will have crossing terms which 



are similar to the case of multichannel Kondo model |27j with the spin index replacing the 
channel index there. The final result is 



U d = (Z^F) 2 ^) 2 (14) 



with the crossing term 



X\{'rii-r2 j ){r li -T 2 



f = Tt7 r XT' Y (15) 

Putting everything together, we find the total contribution of all channels to a given 
term is 

A? 



U T = U U L U c Y[U d = (D ] D i f 1 --Y\D^D i F) 2 ^U*i) 2 . (16) 
i=i 



Following earlier treatments [p6| , p7| , we now consider a generalized resonant level model 



H = H + X[(d+ + 4)^(0) + ^ + (0)(rf T + di)} + ^ + (0)^(0)(n - i), (17) 

where n = (dt + d^)(d^ + d^), with constraint dtd^ + d^d^ < 1, ip + (0), tp(0) a spinless 
fermion. As before, we can expand the evolution operator in terms of A. Again only even 
order terms survive and there are three types of terms (containing only dt, d<\ or dt, or 
mixed. Also, due to the single occupancy constraint any \a > state can be created only from 
the empty state. We can then repeat the same procedure to derive contributions from the 



open line and closed loops. However, there is one important difference, namely, the single 
type of spinless fermion ^(0) is coupled to both and d^, so there are crossing terms 
even for the free propagator (where V = 0). When we switch on the scattering potential at 
moments r ai (the process |0 >^ \a >) ip gets a phase shift 

<5 ' = 2taa- 1 (^) (18) 

and an opposite one at moments r' ai during the process \a |0 >. Summing up the 
contributions from the free propagator, the open line and the closed loops, we find the n-th 
order term of A is given by 

if = (Z^F^-v) 2 . (19) 

If we take 

N r r' 

5 = rr, 2ER 2 = (1--) 2 (20) 
i=i 71 71 

the extended Anderson model H T , given by (D) and (§), and the generalized resonant-level 
model fll~7|) are equivalent to each other via a term-by-term mapping of the perturbation 
expansion. It is expected that they should contain the same low-energy physics. Of course, 
there is an underlying assumption that V is the only scaling parameter for this universality 
class. A similar assumption was made for the single channel Kondo and the two-channel 



Kondo problem with channel asymmetry [27 . 



IV. CORRELATION FUNCTIONS AND PHYSICAL IMPLICATIONS 

Before proceeding we compare first the result of the preceding Section with earlier cal- 
culations using bosonization ]n^,[nj. In the previous work a model very similar to (|1|), (0) 
was considered with additional opposite spin Couloumb interaction Vq c t c cr( n v ~ \) an d 
spin-flip scattering V±J2a c t c ^^d a in the hybridizing channel (see Eq. (6) in |13[). It was 
shown there by using the canonical transformation that the strong coupling Toulouse limit 
is reached for V — > oo, — > 0, V s — \/2N s V s oo, where V s is the Couloumb potential 

10 



in the screening channel and N s is the number of channels. As follows from flTT|), (pT3|), (|i8|) 
and (|20|), we have the same Toulouse limit here. In fact, 5q = tan _1 (7rz/ l / /2) = |, 8' = 
0) 2X^i(^) 2 = 1, meaning that the unitarity limit is reached in both hybridizing and 
screening channels. The earlier calculations correspond to the Born approximation of our 
present result summed to infinite orders. This reconfirms the consistency of bosonization and 
ND approach. However, here we have obtained the mapping to the generalized resonant-level 
model in a broader regime, namely, the unitarity limit should be reached in the hybridizing 
channel (5 = n), but not necessarily in the screening channels (5' ^ 0). Since the Toulouse 
limit is materialized at V ' = 0, the opposite spin scattering is not essential. As seen in 
||13|| , the effect of the V± term is reflected only in the energy difference of a and j3 particles 
(a = -j^(d^ + d±) , P = -^(d^ — di)) which we will include in our following discussion. 

Now we study the physical properties of the generalized resonant-level model (|i7|) 
rewriten as 

H = H + H h + Hi, 

H = Y^ tkipk^k + e a a + a + e^(3 + (3, 

k 

H h = A(a + ^(0)+^ + (0)a), 

Hj = V^ + (0)tlj(0)(a + a-^) (21) 

which should still satisfy the single occupancy constraint a + a + (3 + f3 < 1. First we calculate 
the scattering amplitude 

S{t) =< \e im \ >= exp{C(t)}, 

C(t) =< T{exp[i f dr{H h {r) + i*/(r))]} > c . (22) 
Jo 

There are contributions from the spinless fermion ip and the a particle, both from the closed 
loops. Neglecting the energy level renormalization factors, the long time asymptotics are 
given by 

C(t)~-(^) 2 lnt-(^) 2 lnt, 

7T 7T 

S(t) ~ r^) 2 -^) 2 (l + e-^), (23) 
11 



where 

5 1 = tasr^icpV), 5 2 = tan-^^^). (24) 

with p as the density of states for the spinless fermion ip. Here we have taken the ground 
state energy E = and 82 — f , if e a — 0. There are crossing terms coming from H h and 



Hi of Eq. but they do not contribute to power law singularities in the correlation 

functions. 

Using the asymptotic form for S(t), we can calculate the propagator 

< (3(t)(3 + (0) >~ e-^it)-^ 2 -^ 2 (25) 
and the spin-spin correlation function 

M(t) ~< S*(t)S*(0) >~ cos[(e a - ep)t\(t)-^ )2 -^ 2 (26) 

where S z = \{d^d^ — d^d]) = \(ct + (3 + j3 + a). Assuming Si — 0, S 2 — | , we recover the 
previous result, i.e., M{t) ~ This is a nice and independent check of the correctness 
of calculations in . Moreover, here we have also obtained result away from the unitarity 
limit 5i = 0. Of course, the asymptotic form is valid only for small 81, because the model 



(pT|) cannot be solved exactly. 

As mentioned in the Introduction, the neutron scattering data as well as the static 
measurements in V ' Pd x Cu^- x (x = 1, 1.5) show a power law behavior of the spin-spin corre- 
lation function and the impurity contribution to the specific heat and magnetic susceptibility 



17- 19 1 . Using the conformal invariance for the impurity problem it has been argued that 



the available data are consistent with a critical exponent (T^j A = |, while the value fol- 



lowing from the bosonization calculation [O] was = | which is rather big compared 



with the experimental value. The deviation from the unitarity limit 5i = will reduce this 



value. To fit the data we need to assume 5± = 7ry5/12 which is still rather big to justify 
the applicability of our asymptotic expansion. Nevertheless, the correction is in the right 
direction. 
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V. CONCLUDING REMARKS 



Using the ND and AY perturbation expansion and mapping onto a generalized resonant- 
level model we have reconfirmed and extended earlier results on non-FL behavior in an 
extended Anderson model with additional screening channels. The extension to regions 
away from the unitarity limit of screening electron scattering improves the agreement with 
experiment. However, in view of its physical implications this issue should be further studied 
using other techniques. 

As mentioned in the Introduction, the Anderson model has been also generalized in a 
different way Without going into detailed comparison we briefly comment on the 



study of the Toulouse limit in that model ||29|| . Those authors have correctly pointed out the 
change of sign for the V± term upon bosonization. However, the correlation functions were 
calculated there using the mean field approximation in handling the constraint which missed 
the non-FL exponent. Of course, the physical consequences depend strongly on the positions 
of the renormalized levels e a and ep (see eq. ( p6|) ). If the difference is big, the fast oscillation 
will suppress the power law component in the frequency response. Since the theoretical 
calculation of level renormalization is very difficult, we may count on experimental indication 
which seems to show the existence of remaining degeneracy. The situation here is similar to 
the two channel Kondo case [27|. There the decoupling of d + +d from the conduction electron 
was a signature of non-FL behavior which becomes more apparent in the Majorana fermion 



formulation [31J|. Here the (3 particle decouples, giving rise to a X-ray edge like singularity 
and a residual entropy, as in the case of two-channel Kondo. Of course, this analogy is more 
mathematical than physical. The difference of non-FL behavior in the miulti-channel Kondo 
model and the extended Anderson model in the mixed valence regime, and their possible 
connections have to be further explored. A detailed comparison of the scaling theory with 
numerical Wilson RG studies, as well as a stability analysis of the strong coupling fixed 
point in the extended Anderson remain outstanding issues. 

Finally, we would like to thank G.M. Zhang for an earlier collaboration on this project 
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